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I. INTRODUCTION 

In classical optics, the concepts of the intensity and phase of optical fields have a well- 
defined meaning. That is the electromagnetic field (E) associated with one mode, E = 



O I Aexp{i6), has a well defined amplitude [A) and phase {6). This is not so simple in quantum 



optics where the mean photon number and the phase are represented by non-commuting 
operators and consequently they cannot be defined well simultaneously. In fact, the concept 



H ' of phase is a controversial problem from the earlier days of quantum optics 



m 



2| . In general 

there are three methods of treating this issue [3|. The first one considers the phase as 
a physical quantity in analogy to position or momentum by representing it with a linear 
Hermitian phase operator. The second one involves c-number variables (real or complex) 
in phase spaces or their associated distribution functions, or ensembles of trajectories. The 
third one is the operational phase approach in which the phase information is inferred 
from the experimental data by analogy with the classical analysis of the experiment. Each 
approach has advantage and disadvantage points. Indeed, the interest in the phase properties 
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has been motivated by experimental realization of optical homodyne tomography [^ allowing 
quantum phase mean values to be calculated from the measured field density matrix. 

As is well known squeezed states of light (i.e. the states of light with reduced fluctua- 
tions in one quadrature below the level associated with the vacuum state [51]) have become a 
subject of intensive studies owing to their interesting applications in various devices, e.g. in 
optical communication systems |6|, interferometric techniques |7|, and in an optical waveg- 
uide tap B. Moreover, there are a number of successful experiments producing such hght 
[9[. These states have phase sensitive noise properties and therefore several works have been 



devoted to follow such properties, e.g. for single-mode squeezed states [10|, lUj. On the 
other hand, the superposition principle is at the heart of quantum mechanics. It implies 
that probability densities of observable quantities usually exhibit interference effects instead 
of simply being added. The most significant example reflecting the power of such principle 



are the Schrodinger cat states 1^ , which exhibit various nonclassical properties, even if the 



original states are close to the classical states 



13|]. Based on this principle, a general class 



of quantum states has been introduced as a superposition of displaced and squeezed number 
states [M, ll5| . It is important to point out the basis of this class are the squeezed displaced 



number states 16| which are purely nonclassical states. This class is represented as a single 



mode vibration of the electromagnetic field suddenly squeezed-plus-displaced by a collection 
of two displacements 180° out of phase with respect to each other [ij, [l5|, i.e. 



\r,a,n), = K[b{a) + tb{-a)]S{r)\n), (1) 

where A^ is the normalization constant, D{a) and S{r) are displacement and squeeze operator 
respectively, while a and r are displacement and squeeze parameters; e = |e|e*''^ is a complex 
parameter, and \n) denotes a Fock state. Squeeze and displacement operators are given, 
respectively, by 



5(r) = exp[-(a2-at2)], 



(2) 



D{a) = exp(a^a — aa*), (3) 

where a and a^ are annihilation and creation operators. The normalization constant A^ is 
given by 



|A2|-i = l + |ep + 2|e|exp(-2|tnL„(4|tncos0, (4) 

with 

t = a cosh r + a* sinh r, (5) 

and Ln{-) is the Laguerre polynomiaL The density matrix p of this state can be written as 

p= |AeP(pA/ + P/)- (6) 

The part of the density matrix corresponding to the statistical mixture of two squeezed 
displaced number states is: 

Pm = D{a)S{r)\n){n\S\r)D\a) + \e\^D{-a)S{r)\n){n\S\r)D\-a), (7) 

while the quantum interference part has the form 

Pi = e*b{a)S{r)\n){n\S\r)D\-a) + eD{-a)S{r)\n){n\S\r)D\a). (8) 

This quantum interference part of the density matrix contains information about the quan- 
tum interference between component states D{±a)S{r)\n) and this will be responsible for 
some interesting behaviour of the phase distribution, as we will see. 

For completeness, the physical interpretation of such states can be related to a super- 
position of coherent states formed due to two excitations on particularly excited harmonic 
oscillators 16|. It is clear that these states enable us to obtain generalizations of some 
results given in the literature. It has been shown that these states can be generated, by 
means of the so-called quantum state engineering, and also by means of trapping ions (for 
more details, see refs. l4, llSj). The quantum properties of these states reveal that they can 
exhibit sub-Poissonian statistics, quadrature squeezing and oscillations in photon-number 
distribution. Moreover, the influence of thermal noise on the behaviour of such superpo- 

n 

sition has been considered p/3] showing that the correlation between different oscillators is 

essentially responsible for nonclassical effects similar to Schrodinger cat states. This fact has 

been demonstrated from the behaviour of Wigner function and photon-number distribution. 

In this article we study the phase properties of the superposition ([T]) using Pegg-Barnett 



technique 18j which is most convenient for the current problem. This investigation is orga- 



nized as follows: In section 2 we give the description of the Pegg-Barnett technique and the 



basic relations related to the state ([T]), followed by section 3 where the results are discussed. 
The conclusions are summed in section 4. 

II. BASIC RELATIONS 



Here we give essential background for Pegg-Barnett [ISJ] phase formalism and the basic 
relations for the state under discussion. Their formalism is based on introducing a finite 
(s + 1) -dimensional space \l/ spanned by the number states |0), |1), ..., |s). The physical 
variables (expectation values of Hermitian operators) are evaluated in the finite dimensional 
space ^ and at the final stage the limit s — > oo is taken. A complete orthonormal basis of 
s + 1 states is defined on "^ as 

1 " 

\e^) = ^=^y2expiike^)\k), (9) 

V s + 1 r-" 



V^^.=o 



where 



e„ = eo + ^— , m = o,i,...,s. (10) 

The value of 0o is arbitrary and defines a particular basis of s + 1 mutually orthogonal 
states. The Hermitian phase operator is defined as 

$e = $^e^|e™)(e^|, (ii) 

m=0 

where the subscript shows the dependence on the choice of 60 . The phase states ([9]) are 
eigenstates of the phase operator ( ITTl) with the eigenvalues Gm restricted to lie within a 
phase window between 60 and 27r + Gq. The expectation value of the phase operator ( ITTll 
in a pure state \ip) = Ylm=o^rn\fn), where Cm is the weighting coefficient including the 
normalization constant, is given by 

s 

(^|$,|^) = J]e™|(v^|e™)p. (12) 

The density of phase states is (s + l)/(27r), so the continuum phase distribution as s tends 
to infinity is 



p(e) = iim,^oo^l(e„|V')P 

(13) 
= ^ E^V=oC'mC'% exp[z(m - m')e], 

where Gm has been replaced by the continuous phase variable G. As soon as the phase 
distribution P(G) is known, all the quantum-mechanical phase moments can be obtained as 
a classical integral over G. The phase distribution is normalized such as 

P(G)rfG = 1. (14) 

One of the particular interesting quantities in the description of the phase is the phase 
variance determined by 

{{A^f)=fQ^p{e)de-nQP{e)dQ] . (is) 

As we mentioned earlier the mean photon number and the phase are conjugate quantities 
in this approach and consequently they obey the following uncertainty relation 

((A<l)2)((An)^)>l|([n,<l])p. (16) 

The number-phase commutator appearing on the right hand side of ( IT6|) can be calculated 
for any physical state [l8| as 

([n,<l])=z[l-27rP(Go)]. (17) 



In relation to flT6|) . we can give the notion of the number and phase squeezing 19|,|20|] through 
the relation 

^l(|fl,*l)l 

^ i(A*)l - 1. (19) 

^l{|ft,*l>l 

The values of —1 of these equations means maximum squeezing of the photon number or 
the phase. 
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We shall use above relations to study the phase distribution for the superposition of 
displaced and squeezed number states ([1]). In this case Cm = Cm{r,a,n,e) = {m\r,a,n)^ 
and this quantity can be calculated through the identity 



{m\n, a, r)^ 



+00 



dxT*^{x)Tl^\x,r,a), 



<^) 



where Tn {x,r,a) = {x\n,r,a)^ is the wavefunction of ([T]) having the form 15| 



(20) 



T)^'(x,r,a) 



#^(^)4 <:exp 



f(x^-v^«) 



xH„[e'^(xyf - V2a)] + eexp 



- — [X 



+ V2a)' 



(21) 



H„[e^(x^+V2a)]} 



X 

where a is real, Hn{-) is the Hermite polynomial of order n; u and h are frequency of the 
harmonic oscillator and Planck's constant divided by 27r. Substituting from ( 12T]) together 
with the wavefunction Tm{x) of the Fock state (which can be deduced from (12T]) by simply 
setting r = a = e = 0) into (!20|) and carrying out the integration, using the identity [21 1 



f J_^ dxH„(x)H„(Ax + d) exp(-Mx2 + ex) = exp( 



AM ' 



XI 



M 



M-l\m( I M~A^ \n \^min{m,n) n\m\ 



M 



Emm( 
j=0 



j\(n-jy.{m-jy. 



(22) 



we arrive at |l5 | 



2A 



V/(M-1)(A/-A2) 



^m—j 



TT I cA+2dM 

2y/{M-l)M J ""-^ \2y/M{M-A) J ' 



A / tanh r \ -^r 
/ I \ '^^V 9 ) 



nimicoshr- exp[^(tanhr - 1)] ^-!;;(-'") 



j!(n-j)!(m-j)! I-Vsinh2r-I 



X 



.tanhrl^V^n 



)H. 



=^)[l + (-l)("+-)e; 



(23) 



2 J " " "-■? ^ Vsinh 2r '^ ™~-? ^ Vsinh 2r ' 

where r = aexp(r). Another way for the derivation of fl23p can be found in the appendix 



of 



22]. When r — )■ 0, (123!) reduces to the distribution coefficient for the superposition of 



displaced Fock states as 



C^{r = 0,a,n,e) = A.v/d^exp(-^) Er=o^"^'"^ ^7rSj£5)?^[1 + (-1)"^'"^] 



(24) 



A. 



(-l)"-Pa"+™-2Pexp(-^)[l + (-l)"+'"e]L9-P 



p v« J, 



where i^5^(.) is the associated Laguerre polynomial, and p = niin(n, m) and q = inax{n,m). 
The photon number variance is defined as 



((Ar2)2) = {h^) - {fi)' 



{a)'^a^) + {h) - {hf 



(25) 



where the number operator h = a)a. In terms of (1) the quantities {a)'^a?) and (n) can be 
straightforwardly calculated, however, the explicit forms can be found in |l5|. For complete- 
ness, the phase variance of (1) reads 



((A<|)2) = 2^ + 4Re E„>^. C^{r, «, n, e)C;;,(r, a, n, e 



\-^Y 



[m—m')^ 



(26) 



Rei Em>m' C'm(^' «' n^ e)C;;,(r, a, n, e; 



(-1)" 



■^2 



(m—m') 



The value 7r^/3 is the phase variance for a state with uniformly distributed phase, e.g. the 
vacuum state. Finally, for the future purpose, we give the form of Wigner function of (IT|) 

m 

H/(^)(x,y) = ^^M {|e|2exp[-2(|/2e-2- + e2'-(x + «)2)] 



xLn[4:iy'e-'' + e''{x + af)] + exp[-2{y'e-'' + e^^(x - af)] 



xL„[4(y2e-2r + e^r^^ _ ^)2)] + 2exp[-2(y2e-2^^- + e^^x^)] 



(27) 



X |e|L„[4(|/2e-2'^ + e^^^x^)] cos{4ya - 4>)} . 
Based on the results of the present section we discuss the phase properties of the state 
|1]) in the following section. 



III. DISCUSSION OF THE RESULTS 

For understanding the behaviour of the phase distribution of the states ([1]) it is more 
convenient firstly to study such properties for two subsidiary states which are the superposi- 



tion of displaced number states [23|, and squeezed and displaced number states [16]. These 
two states are of interest because it has been shown for the former states that they can 
exhibit strong sub-Poissonian character as well as quadrature squeezing. The bunching and 



antibunching properties for the latter states have been discussed in 2J|. Investigation of 
such properties for these two kinds of states most clearly illustrates the role of the different 
parameters in the state ([H) in the behaviour of the distribution. 

In the following two subsections we discuss the phase distribution determined by ( TT3l) . 
and the phase variance and phase squeezing, respectively. For simplicity we restrict our 
investigation to real values of a and e. 

A. Phase probability distribution 

In general we find that there are two regimes controlling the behaviour of the phase 
distribution for the states under discussion which are a >> 1 and a < 1 provided that 
n is finite. For the first case, the even- and odd-cases (i.e. e = 1 and —1) provide similar 
behaviours and there is one-to-one correspondence between the components of density matrix 
([6]) and that corresponding to the curves. In other words, pj\/ is responsible for the central 
behaviour (i.e. around O = 0) of the phase distribution and pi is responsible for the lateral 
behaviour, i.e. when — t- ±7r (see Fig. 2a). Nevertheless, all these features are washed 
out in the second regime and the behaviour dramatically change. Such behaviours can 
be understood well if we turn our attention to the behaviour of the PF-function where the 
interference in phase space may be seen clearly. More illustratively, W^-function ( 1271) includes 
three terms, the first two terms representing the statistical mixture of squeezed displaced 
number states and the third one is the interference part. When a < 1 (as shown in Fig. la 
and b for superposition of displaced number states) the contributions of these components 
are comparable so that even- and odd-cases distributions have different shapes. On the 
other hand, when a >> 1 the structures of the VT-function for even- and odd-cases are 
almost similar, i.e. they include two symmetrical peaks originated at ±a and interference 
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fringes are in between. Of course, such situation is still valid if squeezing in the displaced 
superimposed number state in optical cavity is considered, however, the peaks will be then 
stretched. This behaviour is reflected in the behaviour of the phase distribution, as we will 
see. It should be stressed that the greatest degree of the nonclassical behaviour for the 
superposition (1) actually occurs for rather low values of a |lJ, Il5| . 

(la) 

(lb) 




(a) (b) 

FIG. 1: ly-function against x and y for n = 1, a = 1 and for a) e = 1; b) e = —1. 

We start our discussion by focusing the attention on the behaviour of the superposition of 
displaced number state. It is important to mention that the properties of displaced number 
states D{(x)\n) have been given in 25| and their phase properties in 26|. In fact the phase 
properties of such states are interesting since they connect the number state, which has no 
phase information, and coherent states, which play the boundary role between the classi- 
cal and nonclassical states and always exhibit single peak structure of phase distribution. 
The phase distribution of displaced number states exhibits nonclassical oscillations with the 
number of peaks which are equal to n + 1. This result is interpreted in terms of the area 
of overlap in phase space 261]. Unfortunately, the situation is completely different for the 
superposition of displaced number state (see Fig. 2a and b for shown values of the param- 
eters). From Fig. 2a it is clear that the quantum interference between component states 
D[(y)\n) and D{^—a)\n) leads to the lateral nonclassical oscillations which become more pro- 
nounced and narrower as n increases. However, the statistical mixture of displaced number 
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states shows n + 1 peaks around 9 = 0. We have not presented the phase distribution for 
the odd-case since it is similar to that of the even-case. Comparing this behoviour with 
that of even (or odd) coherent states, we conclude that the number states in the superposi- 
tion make the phase information more significant [Ij. Turning our attention to the Fig. 2b 
where a = 1, we can see that the distribution of displaced number states exhibits two-peak 
structure as expected for n = 1 [261, however, the distribution of even-case displays one 
central peak at O = and two wings as O — ?■ ivr, and finally the distribution of odd-case 
provides four peaks. That is the distribution here is irregular, however, more smoothed than 
before and the structure of the statistical mixture (central part) is modified by the action 
of the interference term arising from the superposition of the states; this was the case for 
VT-function. 

(2a) (2b) 

0.50n 0.40n 



0.30- 
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^ - 



0. 10- 
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4. 00 
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-4.00 



.00 0.00 2.33 

Theta 



4. 00 



(a) (b) 

FIG. 2: Phase distribution P{Q) for superposition of displaced number states and for a) a = 
2,e = 1 , n = 1 (dashed curve), and 2 (solid curve); b) a = l,n = 1 and e = (bell-centered 
curve), 1 (solid curve), and —1 (dashed curve). 



Before discussing the phase properties of the displaced and squeezed number states it is 
reasonable to remind the behaviour of the well known squeezed states. As known for squeezed 
states with non-zero displacement coherent amplitude, the phase distribution exhibits the 
bifurcation phenomenon. In this phenomenon the single peak structure of the coherent 
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component is evolved into two peaks structure with respect to both a (for large fixed value 
of squeezing parameter r) and r (when a takes fixed value) ll|]. This phenomenon has 
been recognized as a result of the competition between the two peaks structure of the 
squeezed vacuum state and the single-peak structure of the coherent state. For squeezed 
and displaced number states such phenomenon cannot occur due to the effect of the Fock 
state which replaces the initial peak (r = 0) for coherent state by a multi-peak structure, 
i.e. by ra + 1 peaks (see Fig. 3a for shown values of the parameters). From this figure 
one can observe that there is a three-peak structure corresponding to the case n = 2 of 
displaced Fock state. The height of the central peak (i.e. at G = 0) is almost the same and 
equals (l/27r)| X]m=o^™l^- That is the central value of the phase distribution P(6 = 0) is 
insensitive to squeezing provided that r is finite. However, the lateral peaks undergo phase 
squeezing as r increases, i.e. the peaks become narrower. 

Now we can investigate the behaviour of the superposition of displaced and squeezed 
number states (see Figs. 3b-d for shown parameters). Figs. 3b and 3c are given for the 
first regime for even- and odd-cases, respectively. From Fig. 3b we can see that the initial 
oscillations are increased compared with Fig. 3a as a result of the interference in phase 
space. Further, the initial lateral peaks can evolve in the course of increasing r to provide 
bifurcation shape, i.e. the distribution curve undergoes a transition from single- to a double- 
peaked form with increasing r; however, the central peak is almost unchanged with increasing 
r. This peak splitting is connected with the squeezed states HI]. Further, as the number of 



peaks increases for r > 0, the distribution becomes more and more narrower. On the other 
hand, the phase distribution of the odd-case is quite different as we have shown before where 
the initial peaks are not significantly changing as in the even-case. More precisely, initial 
distribution becomes broader for a while and suddenly (at r ~ 0.5) breaks off to start to 
be a narrower distribution for later r. It is clear that in this regime the state ([1]) becomes 
more and more nonclassical [15| and the even- and odd-cases are distinguishable. Further, 
the phase distribution of the even-case is more sensitive with respect to squeezing than the 
odd-case. Nevertheless, for the second regime a >> 1 we noted that the phase distribution 
carries at least the same initial information regardless of the value of r (see Fig. 3d for the 
shown values of the parameters). 
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FIG. 3: 

(e,a,n) 
tively. 



(c) (d) 

Phase distribution P{Q) for superposition of displaced and squeezed number states for 
= (0, 1, 2), (1, 1, 2), (— 1, 1,2) and (1,2,2) corresponding to the cases a,b,c and d, respec- 



B. Variance, amplitude and phase squeezing 



In this subsection we investigate the behaviour of the phase variance, and amplitude and 
phase fluctuations following f l26|) . flTSj) and f lT9|) . respectively. We start our discussion by 
analysing the behaviour of the superposition of displaced number states. For this purpose 
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Figs. 4a and 4b are shown for the phase variance, and amphtude and phase fluctuations, 
respectively. It is seen that in general the phase variance starts from the value 7r^/3 (the 
vacuum state value) and returns back to it when a is large, but through different routes. To 
be more speciflc, the phase variance of displaced number states starts from the value for the 
vacuum, goes to a minimum, and then comes again to vr^/S. However, the behaviour of the 
superposition of displaced number states takes different ways to arrive at the same result, 
i.e. it starts from 7r^/3 as before, goes to the maximum value and eventually comes back 
to the value of vacuum. The comparison of the two cases shows the role of the quantum 
mechanical interference between state components. Further, as n increases, the oscillations 
in the variance become more pronounced. Comparison of the behaviour of even- and odd- 
cases shows that they are different only over the initial short interval of a, i.e. when a is 
small, and this agrees with what we have discussed earlier. So we can conclude that for 
intensities high enough of the coherent fleld, the variance of the phase is approximately 
randomized. The route to this randomization is dependent on the choice of e. With respect 
to the amplitude and phase fluctuations, we can note from (ITTl) that these quantities depend 
not only on the intensity of the fleld, but also on the choice of the reference angle Bq. We 
have chosen here 0o = — vr, where the mean value |([A^, ^d\)\ of squeezed displaced number 
state approaches unity. Fig. 4b has been obtained to illustrate the parameters S^ and Sq 
which provide information about the degree of squeezing in h and $9. One can observe from 
this flgure that when e = (displaced number state) and a — )■ 0, the parameter Sj^ tends to 
—1, which means that the number state is 100% squeezed with respect to the operator h. 
This situation is expected since ((An)^) = for the number state. Further, the larger the 
number of quanta is the shorter is the interval over which Sfq is squeezed. Also when a >> 1, 
squeezing in Sq is remarkable, whereas S^ becomes unsqueezed. This result can be deduced 
from the behaviour of the phase variance (see Fig. 4a), where ((A<l>)^) ^ at a ~ 4 and 
this should be connected with maximum squeezing in Sq at this point. Such behaviour of 
Sfq and Sq conflrms the fact that the number of photons and phase are conjugate quantities 
in this approach. As is known displaced number states are not minimum uncertainty states 
and the variances for the quadrature operators never go below the standard quantum limit. 
Moreover, they may exhibit sub-Poissonian statistics for the range a^ < 1/2 [13^. However, 
there is no relation between the sub-Poissonian statistics and the fluctuation in the amplitude 
or the phase. This fact has been shown before for the down-conversion process with quantum 



14 
pump where the signal mode can exhibit amphtude squeezing and at the same time it is 



super- Poissonian [27 1. 



We proceed by discussing the behaviour of the superposition states (long-dashed and 
circle-centered curves) in Fig. 4b where the interference in phase space starts to play a role. 
We noted (from our numerical analysis) that only the even case can provide squeezing in Sn 
with maximum value at the origin and squeezing interval larger than that discussed before. 
Indeed, this maximum value is related also to that of number state at a = 0. It should be 
stressed here that a = for the odd-case may lead to a singularity. As we mentioned before 
the superposition of displaced number states can exhibit strong sub- Poissonian character as 
well as quadrature squeezing 23|. Now we can illustrate the role of the squeezing in the 
superimposed displaced number states optical cavity with respect to variance, amplitude 
and phase fluctuations. It is obvious, when squeezing is considered, that the initial value (at 
a = 0) of the phase variance is shifted since we have initially squeezed number state which 
is providing phase information. However, when a is large and r is finite or also r is large, it 
can be proved simply that the coefficient Cmir, a, n, e) vanishes and consequently the phase 
variance tends to vr^/S (becomes randomized). Moreover, the routs are here similar to those 
of Fig. 4a. On the other hand, squeezing could be seen only in Sq for e = (see squared- 
centered curve in Fig. 4b). Furthermore, comparison of the short-bell-centered curve (of 
displaced Fock state) and squared-centered curve reveals that squeezing parameter reduces 
the amount of squeezing in Se, too. This means that the superposition of displaced and 
squeezed number states which provide quadratures squeezing has less information about 
amplitude and phase fluctuations. 

IV. CONCLUSIONS 

We have discussed the phase properties of the superposition of squeezed and displaced 
number states from the point of view of the Pegg-Barnett Hermitian phase formalism. More- 
over, the phase variance, the photon number and phase fluctuations have been discussed, 
too. The results have been explored graphically. 

We have shown that there are two regimes controlling the behaviour of the phase for 
such superposition depending on whether the superposition is macroscopic (a >> 1) or 
microscopic (a < 1). This fact has been illuminated by means of the Wigner function. In 
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FIG. 4: a) Phase variance of a superposition of displaced number states against a for (e, n) = 
(0,3) (long-dashed curve), (1, 1) (bell-centered curve), (1,3) (solid curve) and (—1,3) (short-dashed 
curve). The solid straight line is corresponding to the phase variance of vacuum, b) Amplitude 
and phase fluctuations of a superposition of displaced and squeezed number states against a for 
{t^n^r^ = (0,1,0) ( Sn solid curve and Sq short-bell-centered curve), (0,2,0) ( Sn short-dashed 
curve and Sq long-bell-centered curve), (1, 1, 0) {Sn long-dashed curve), (1, 2, 0) {Sn circle-centered 
curve) and (1, 2, 0.5) {Sq squared-centered curve). The solid straight line is the bound of squeezing. 

the first regime, the even- and odd-cases (i.e. e = 1 and —1) give similar behaviours and the 
structure of the density matrix is remarkable in figures. All these facts are washed out in 
the second regime where the behaviour becomes irregular, however, smooth. In general we 
noted that the higher the number of quanta is, the more peaks the distribution possesses. 
Influence of the squeezing parameter could be recognized in the second regime where the 
distribution exhibits peak-splitting and peak-narrowing and this is in contrast with the fist 
regime. 

For the phase variance we conclude that it asymptotically goes to the value 7r^/3 of the 
uniform distribution when either a or r is large, but through different routes. We have shown 
also that this superposition can exhibit photon number fluctuations and phase fluctuations. 
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